

// z0 vs hitY
etf->Draw("sqrt(mu0*exp(0.0018*(dcTaU1_wireHitPos[1]+30.))/(0.004*pow(atof1-116., 2.)+5.))*53.9046-0.480279:dcTaU1_wireHitPos[1]>>h0(36, -45., 45, 1000, 7., 43.)", "dcTaU1_wireHitPos[1]!=-999.&&atof1>0.", "colz")

// aoz vs bg0
etf->Draw("1.6058799/abg0*(1.+0.000194047*dcTaU1_wireHitPos[0]-0.0011855*xf1)+0.19894667:atof1>>h5(1000, 120, 135., 1000, 1.6., 2.27)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&az0>22.&&az0<23.", "colz")
etf->Draw("1.6058799/abg0*(1.+0.000194047*dcTaU1_wireHitPos[0]-0.0011855*xf1)+0.19894667:atof1>>h5(500, 120, 135., 500, 1.6., 2.27)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&az0>27.5&&az0<28.5", "colz")

// aoz vs tof1
etf->Draw("1.6058799/ETFMath::BetaGamma(86.329723/(atof1))*(1.+0.000194047*dcTaU1_wireHitPos[0]-0.0011855*xf1)+0.19894667:atof1>>h5(1000, 120, 135., 1000, 1.6., 2.27)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&az0>22.&&az0<23.", "colz")

// all raw
 etf->Draw("sqrt(mu0*exp(0.0018*(dcTaU1_wireHitPos[1]+30.))/(0.004*pow(atof1-116., 2.)+5.))*53.9046-0.480279:1.6058799/ETFMath::BetaGamma(86.329723/(atof1+10))*(1.+0.000194047*dcTaU1_wireHitPos[0]-0.0011855*xf1)+0.19894667>>h6(1000, 1.9, 2.4, 1000, 0., 50.)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&atof1>0.&&86.329723/(atof1+10)<1.", "colz")
// original
etf->Draw("az0:aaoz0>>h2(1000, 2.0, 2.4, 1000, 10., 41.)", "", "colz");

// final - 20240609_1950.dat_0.root
etf->Draw("sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*53.9046-0.480279:1.6058799/ETFMath::BetaGamma(86.33/(atof1+13))*(1.+0.0003*dcTaU1_wireHitPos[0]-0.00105*xf1-0.000002*xf1*xf1)+0.19894667>>h7(1000, 1.9, 2.5, 1000, 0., 50.)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&86.33/(atof1+13)<1.&&atof1>0.", "colz")

etf->Draw("sqrt(mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.))/2000.)/(0.004*pow(atof1-116., 2.)+5.))*53.9046-0.480279:dcTaU1_wireHitPos[1]>>h0(36, -45., 45, 1000, 14., 30.)", "dcTaU1_wireHitPos[1]!=-999.&&atof1>0.", "colz")

// z0 vs hitY
etf->Draw("sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103:dcTaU1_wireHitPos[0]>>h7(36, -45., 45., 1000, 5., 40.)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&86.33/(atof1+13)<1.&&atof1>0.&&sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103>0.", "colz")
// z0 vs x_f1
etf->Draw("sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103:xf1>>h8(52, -52., 52., 1000, 5., 40.)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&86.33/(atof1+13)<1.&&atof1>0.&&sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103>0.", "colz")
// z0 vs tof1
etf->Draw("sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103:atof1>>h9(1000, 115., 135., 1000, 5., 40.)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&86.33/(atof1+13)<1.&&atof1>0.&&sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103>0.", "colz")

// run: 20240609_1950.dat_0.root
// calibrated - both aoz and z - 23:03, Jul 4th, 2024
etf->Draw("sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103:1.6058799/ETFMath::BetaGamma(86.33/(atof1+13))*(1.+0.0003*dcTaU1_wireHitPos[0]-0.00105*xf1-0.000002*xf1*xf1)*1.25631-0.249302>>h7(1000, 1.9, 2.6, 1000, 5., 40.)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999.&&dcTaU1_wireHitPos[1]!=-999.&&86.33/(atof1+13)<1.&&atof1>0.", "colz")

// xf1 vs xETF -- per Z
etf->Draw("xf1:dcTaU1_wireHitPos[0]>>h0(36, -45., 45., 52, -52., 52.)", "xf1!=-999.&&dcTaU1_wireHitPos[0]!=-999. && sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103>12.-0.5      && sqrt((mu0v[2]*exp(0.0015*(dcTaU1_wireHitPos[1]+30.))+mu0v[3]*exp(0.0022*(dcTaU1_wireHitPos[1]+30.)))/2000./(0.0016*pow(atof1-100., 2.)+5.))*50.7183-1.34103<12+0.5", "colz")

